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Abstract 

We present an overview of recent works on the statistical description of turbulent 
flows in terms of probability density functions (PDFs) in the framework of the 
Lundgren-Monin-Novikov (LMN) hierarchy. Within this framework, evolution 
equations for the PDFs are derived from the basic equations of fluid motion. 
The closure problem arises either in terms of a coupling to multi-point PDFs or 
in terms of conditional averages entering the evolution equations as unknown 
functions. We mainly focus on the latter case and use data from direct numerical 
simulations (DNS) to specify the unclosed terms. Apart from giving an introduc- 
tion into the basic analytical techniques, applications to two-dimensional vorticity 
statistics, to the single-point velocity and vorticity statistics of three-dimensional 
turbulence, to the temperature statistics of Rayleigh-Benard convection and to 
Burgers turbulence are discussed. 

Nous presentons un apergu des travaux recents sur la description statistique des 
ecoulements turbulents en terme de fonctions densite de probability (PDFs) dans 
le cadre de la hierarchie de Lundgren-Monin-Novikov (LMN). Dans ce cadre, 
des equations devolution pour les PDFs sont derivees a partir des equations 
fondamentales decrivant la dynamique des fluides. Le probleme de fermeture se 
pose soit sous la forme d'une couplage aux PDFs multipoints, soit sous la forme 
de moyennes conditionnelles qui entrent dans les equations devolution en tant 
que fonctions inconnues. Nous nous concentrons principalement sur le dernier 
cas et utilisons les donnees des simulations numeriques directes (DNS) pour 
calculer les termes non fermes. Nous donnerons done une introduction de base 
aux techniques analytiques. Ensuite, nous presenterons quelques applications, en 
particulier a la statistique de la vorticite en deux dimensions, aux statistiques 
de vitesse et de vorticite en un point en trois dimensions, a la statistique de 
la temperature dans la convection de Rayleigh-Benard et a la turbulence de 
Burgers. 



Email address: mwilczek@uni-muenster.de (M. Wilczek 1 ) 



Preprint submitted to Elsevier 



October 23, 2012 



Contents 



1 Introduction [3] 

1 . 1 Outline of the Review [5] 

2 Kinetic Equations: Concepts and Methods [6] 

2.1 Probability Distributions and Random Fields [5] 

2.2 Statistical Symmetries [5] 

2.3 LMN Hierarchy for the Vorticity in Two-Dimensional Turbulence QT| 

2.4 The Method of Characteristics QJ] 

2.5 From the Eulerian to the Lagrangian Frame [TB] 

3 Inverse Cascade in Two-dimensional Turbulence 1171 

3.1 Two-Point Kinetic Equation [T7] 

3.2 Gaussian Approximation to the Conditional Velocity Field .... [18] 

3.3 Failure of the Gaussian Approximation Q]|] 

4 Velocity Statistics in Three-Dimensional Turbulence 1211 

4.1 LMN Hierarchy for the Velocity |2T] 

4.2 Single-Point Velocity Statistics \M 

5 Vorticity Statistics in Three-Dimensional Turbulence 1271 

5.1 LMN Hierarchy for the Vorticity [27] 

5.2 Single-Point Vorticity Statistics [Ml 

6 Application to Turbulent Rayleigh-Benard Convection 1321 

6.1 LMN Hierarchy for Temperature and Velocity [33J 

6.2 Numerical Results [35] 

7 Burgers Turbulence 1361 

8 Conclusions 1381 



2 



1. Introduction 

The challenge in fundamental turbulence research is the prediction of char- 
acteristic flow features starting from the basic equations of fluid motion. For 
flows in the turbulent regime predictive results are rather sparse and, despite of 
important applications in science and engineering, basic properties of turbulent 
flows remain poorly understood starting from first principles pQ. 

The same applies to statistical hydrodynamics, which aims at a description 
of turbulent flows based on a statistical treatment of the basic hydrodynamic 
equations [2J|5]. In particular, solving the statistical equations derived from 
the basic hydrodynamic equations has only been successful in but a few cases 
with Kohnogorov's four- fifth law being one of the most prominent examples (4]. 
The main reason for this is that the Navier-Stokes equation, which describes the 
fluid motion mathematically, is both nonlinear and nonlocal. As a consequence, 
statistical relations derived from this equation are in general unclosed, such 
that a statistical description without further ad hoc assumptions usually leads 
to an infinite set of statistical equations or, in a more condensed notation, to 
functional equations [2J [3] . The investigation of these statistical equations and 
the development of physically reasonable closures may be seen as one of the 
central issues of theoretical turbulence research. 

The statistical investigation of turbulent flows has a long history, with its 
origin in the works of Reynolds [5] , who formulated moment equations for the 
mean and fluctuating velocity. Already this pioneering work has revealed the 
fundamental problem of statistical turbulence research, the closure problem, 
showing up in terms of the turbulent Reynolds stresses. The infinite hierarchy 
of moment equations was first formulated by Friedmann and Keller [5] and has 
become the major line of research during the first half of the twentieth century. 
The challenge consists in adding further exact relations to Kohnogorov's four- 
fifth law, which allows one to formulate exact or approximate closures for the 
hierarchy 015]. The formulation of an ad hoc closure for a higher-order moment 
equation eventually may lead to a realizability problem, e.g. the generation of 
negative energy. This realizability problem arises in the naive Gaussian closure 
formulated by Millionshchikov [9]. 

Instead of formulating ad hoc closures of the Friedmann-Keller hierarchy, 
it is possible to study the statistics in the framework of the Lundgren-Monin- 
Novikov hierarchy, which often also is referred to as "PDF methods" or "kinetic 
equations". Work on this has been initiated in the sixties of the last century 
by Lundgren, Monin and Novikov |10[ [TTJ [12]. In this theoretical framework, 
evolution equations for the probability density functions of turbulent observables 
are obtained from the equations of motion of the fluid. Generally speaking, 
kinetic equations are conservation laws for the probability distributions of the 
vorticity or velocity field at a number of spatial points and time instants. They 
are obtained by relating the dynamics of the TV-point probability distributions 
to the dynamics of the underlying classical field theory. The result is an infinite 
hierarchy of evolution equations for the TV-point distribution functions of the 
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form 



0_ 
di 



}n — Lw/at + Ljv,jv+i/jv+i ■ 



(1.1) 



Here, /jv and Jn+i are the iV and (N + l)-point PDFs, which are related by 
operators Ljy and Ltv^+i, where the latter explicitly highlights the coupling 
to the higher-order probability density. While Lundgren and Monin focused 
on the velocity statistics, Novikov studied the statistics of the vorticity field. 
Of course, unclosed terms arise also in this framework, which can be expressed 
either in terms of the statistics of an increased number of points in space or in 
terms of conditional averages, which then enter the PDF equations as unknown 
functions. Lundgren and Monin mainly followed the former way and investigated 
the arising hierarchy with respect to possible analytical closures. Additionally 
Lundgren pointed out interesting analogies to the BBGKY hierarchy of kinetic 
gas theory [T3| [2]. Novikov emphasized the possibility to introduce conditional 
expectations. A particularly interesting paper is the one of Ulinich and Lyubimov 
|15j . which shortly appeared after the works of Lundgren, Monin and Novikov. 

It can be shown that the LMN approach is completely equivalent to the 
statistical description of turbulence by moment equations. In fact, the Friedmann- 
Keller equations for the A-point-single-time correlation functions of the field 
can be obtained by deriving the evolution equations for the expectations 

(u;(x 1 ,t)...u>(x N ,t)) = Jduji...duj N uji...u} N f({ui,xi},t) (1.2) 

in a straightforward manner, here exemplified for the pseudo-scalar vorticity 
of two-dimensional turbulence. The infinite set of evolution equations for the 
correlation functions, which frequently are also denoted as the Hopf equations, 
can be summarized in a condensed notion in the Hopf-functional equation for 
the characteristic functional 



In this sense, the LMN hierarchy is equivalent to Hopf 's functional equation. 
In fact, Monin showed that the hierarchy of PDF equations corresponds to the 
Hopf functional equation formulated at discrete points in space. 

The motivation for the study of the LMN hierarchy is similar to the one for 
the examination of the Friedmann-Keller equation. The aim is to find a closure of 
the chain of equations for low values of N, which allows one to assess the statistics 
of the field restricted to N points in space. The advantage of formulating closures 
for the LMN hierarchy compared to the Friedmann-Keller hierarchy consists 
in the fact that the formulation of realizable closures, i. e. closures which do 
not generate spurious results like negative energy spectra, is more transparent 
compared to closures for the Friedmann-Keller hierarchy. This has been stressed 
essentially by Lundgren and, later on, has been advocated by Pope in his work 
on turbulent combustion [16| . Furthermore, there is a (formal) analogy to the 
kinetic theory of many-particle systems and the derivation of kinetic equations 




(1.3) 
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using the hierarchy of Bogoljubov, Born, Green, Kirkwood, Yvon (BBGKY 
hierarchy O E]), and it is tempting to apply decoupling procedures developed 
in the kinetic theory of gases and plasmas to the LMN hierarchy. 

Besides the possibility to formulate realizable closures, the LMN hierarchy 
has another appealing property, which is related to progress in the held of direct 
numerical simulations (DNS). This line of research actually has been initiated by 
Novikov who started his work with the infinite set of evolution equations while 
later focusing on the conditional vorticity field combining analytical closures and 
numerical investigations of the field. The introduction of conditional probability 
distributions 

Jn+i = Pn+i,nIn (1-4) 

allows one to formulate (formally closed) kinetic equations for the iV-point 
probability distribution: 

d 

w,Jn — [Lat + M N ]f N , Mjy = Ljv ,n+iPn+i,n (1-5) 
ot 

Progress in DNS allows one to assess the conditional PDFs for the low N 
case, especially for N = 1,2. An interesting conclusion from Novikov 's work is 
that the statistical equations describing the vorticity are incompatible with the 
assumption that the vorticity field can be represented by a Gaussian random 
field [12,. This shows that a statistical description of this non-equilibrium system 
has to rely on the usage of non-Gaussian random fields. This is certainly one of 
the main conceptual challenges of modern turbulence research. 

PDF equations in combination with conditional averages have also been 
established by Pope for the investigation of reactive flows and combustion, based 
on the observation that the reaction kinetics is local in space, and, hence does 
not require the formulation of a closure in the corresponding PDF equation. 
Combined with stochastic modeling, a whole branch of research with engineering 
applications has emerged. Although the starting point of Pope's investigations 
is very close to the work of Lundgren, his line of research aims at modeling and 
applications [TBI H], especially with respect to single-point closures. 

Theoretical research on the LMN hierarchy remains an active area of research 
up to today. Recent work by Yakhot deals with analytical closures and the 
application to scalar fields and convection |171 fTHl fT5] . Work on the statistics 
of velocity fields has been pursued by Hosokawa as well as Tatsumi et al. , with 
the aim of formulating analytical closures of the PDF equation (2UJ ED 121] • An 
especially interesting step here was taken by Tatsumi with his so-called cross- 
independence hypothesis, which neglects statistical correlations of the increments 
and sums of velocities at two points in space, yielding a non-trivial closure 
to the statistical equations [23 . Boffetta et al. [24J used DNS to formulate a 
closure for the statistics of velocity increments of two-dimensional turbulence. 
Ching [25J considered applications to convection, and furthermore, together with 
Pope, formulated equations for the PDFs containing conditional averages [261 12"7] . 
These relations are based on statistical symmetries, and are reviewed below. 

Finally, we mention the work of Polyakov [28 on the statistical formulation 
of the forced Burgers turbulence. He formulated the evolution equation for the 
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single-point PDF and a closure of the dissipative term. His work stimulated 
a variety of theoretical investigations, summarized in the review of Bee and 
Khanin [29] . 

Whereas the approaches based on the formulation of kinetic equations aim at 
assessing the statistics of the turbulent helds starting from exact relations derived 
from the basic fluid dynamics equation, a completely different approach has been 
pursued in the works of Friedrich, Peinke and coworkers |301 15T1 1521 1551 151] . In 
this approach a phenomenological theory has been developed, which is based on a 
Fokker-Planck equation for the PDFs of the velocity increments. The evolution of 
the velocity increments v(r, t) = u(x + r,t) — u(x, t) is considered as a Langevin 
process in scale r. Empirical investigations of this process prove Markovian 
properties on scales larger than the Einstein-Markov scale [35] (which is of 
the order of the Taylor length). Furthermore, estimation of the corresponding 
Kramers-Moyal coefficients allows one to formulate a Fokker-Planck equation for 
the transition PDFs as a function of scale. In this way one gets access to the 
iV-scale velocity increment PDFs and obtains a phenomenological description 
of the phenomenon of intermittency. One may hope that eventually the kinetic 
equation approach to turbulence based on the LMN hierarchy, which is the 
topic of the present review, may lead to a derivation of the phenomenological 
Fokker-Planck theory of intermittency of the turbulent cascade. 

1.1. Outline of the Review 

The aim of the present review is to summarize a number of recent works 
applying the statistical equations of the LMN hierarchy to some basic turbulent 
flows. It is structured as follows. In the next section we summarize basic concepts 
and methods of the kinetic equation approach to turbulent fields. Then, we 
review recent work on the two-point statistics of two-dimensional turbulence, 
and single-point statistics of velocity and vorticity of fully developed three- 
dimensional turbulence. As an example of a confined turbulent flow we discuss 
Rayleigh-Benard convection with emphasis on the boundary layer statistics of 
the turbulent field. We end up with some analytical considerations of Burgers 
turbulence before we conclude. 

2. Kinetic Equations: Concepts and Methods 

In this subsection we introduce the basic definitions and some relations that 
are necessary for the statistical description of classical random fields. For the 
example of the vorticity statistics of two-dimensional turbulence we show how 
these techniques can be utilized for the derivation of the evolution equations 
of the LMN hierarchy. For a very readable account on this topic we refer the 
reader to the paper by Lundgren [TU] or the textbook by Pope [I]. 

2.1. Probability Distributions and Random Fields 

The starting point is the definition of the so-called fine-grained probability 
distribution for the vorticity field at a base point X\ at time t, 

fi(u)i,xi,t) = S(wi - u}(xi,t)) . (2.1) 



G 



Here u(xi, t) denotes a realization of the vorticity field and u>i is the correspond- 
ing sample space variable. Therefore, f\ is a PDF with respect to oji and a 
function with respect to the variables X\ and t. The fine-grained PDFs, loosely 
speaking, provide a probabilistic description of a single realization of the random 
field. In order to obtain the full PDFs from the fine-grained PDFs, one has to 
perform an ensemble average over all possible realizations of the flow leading to 

/i(wi,xi,t) - (/i(wi,xi,t)) = (Sfa-ufrut))), (2.2) 

where the brackets denote the averaging process. Both definitions can be easily 
extended to N points yielding 

f N (uJl,Xl, . . . ,UJ N ,X N ,t) = (f N (ui,Xl, . . . ,LO N ,X N ,t)) 

= (j[6(to i -Lj(x i ,t))y (2.3) 

In general the PDFs have to fulfill some constraints like the normalization 
property 

JdU! . . . JdLJ N f N = 1 (2.4) 

and the reduction property 



/ 



du^ f N = f N _! . (2.5) 

Due to the fact that we investigate physically meaningful fields, the PDFs are 
also subject to additional constraints like the separation property 

lim f 2 (uJi,x 1 ,uj2,X2,t) = fi(uii,x 1 ,t)f 1 (u} 2 ,x 2 ,t) (2.6) 

which means the vorticity at two different points should become statistically 
independent when the distance between them is large enough, and the coincidence 
property 

lim f 2 {u 1 ,x 1 ,LJ 2 ,x 1 ,t) = 6(wi - w 2 )/i(wi,xi,i) , (2.7) 

|as 2 -»i|->-0 

which describes the fact that if the separation between two points becomes 
infinitesimally small, the vorticities uji and cj 2 should also converge to the same 
value, i.e. the vorticity field is smooth on small length scales. 
Beside the joint PDFs we can define the conditional PDF 

/ I .n fN(ul,X 1 ,U> 2 ,X 2 ,...,UJ N ,X N ,t) 

P(uj 1 ,x 1 \uj 2 ,x 2 , . . .,uj N ,x N ,t) = — — (2.8) 

/JV-l(W2, x 2 , . . . ,0J N , X N ,t) 

which is the probability of finding u)\ at X\ for given lo 2 , . . . , ujn at x 2 , . . . ,Xn- 
Now we can also define the conditional averages 

(u}(xi,t)\w 2 ,X2,. ■ ■ ,uj N ,x N ,t) = J Aunoj\p{uJi,xi\w 2 ,x 2 ,...,w N ,x N ,t) 

(2.9) 
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which allows us to write down relations such as (see Appendix H in j3] for details) 



= (aj(xi,t)\u> 2 ,x 2 ,-. ■ ,0J N ,x Nl t)f N -i(uj2,X2, . . .,w N ,x N ,t) (2.10) 

that will play a crucial role for the derivation of the LMN hierarchy. 

To keep the notation as compact as possible we will from now on skip the 
indices of the PDFs, because the order of the PDF becomes clear by its arguments. 
Additionally we introduce the abbreviation 



f({ui,xi},t) := f N (u)i,Xi, . . .,u N ,x N ,t) 



(2.11) 



for the A-point PDF where {uii, xi} denotes the set of N points u\, x\, . . . , ojn, xn- 
In the following we will also frequently use the notation /(a/, x' , {uji, x{\ : t) for 
the A + 1-point PDF. 



Since conditional averages like e.g. in eq. (2.101 are central quantities, it 



would be helpful to derive some general analytical expressions for them. In case 
the characteristic function of the N + 1-point distribution f(to', x', {w;, xi}, t) is 
known the conditional averages can be expressed by 



(u',x'\{u h xi})f({LJi,xi},t) 
1 d 
i da' 



7 W[ a',x', 



1 

i duji 



,x, 



f({oj h xi},t) (2.12) 



a'=0 



where W(a' , x' , {a>i, xi}) is the negative logarithm of the characteristic function 
known as cumulant generating function. For a multivariate Gaussian PDF the 
conditional expectation can be expressed in terms of the correlation functions. 
The result is 



(u/, x'\{uji, xi}) = C(x! - x k )C 1 (x k - x m )uj n 



k.m 



(2.13) 



where C(x! — x k ) stands for the two-point vorticity correlation functions, and 
C~ 1 (xi — Xj) denotes the (i,j)-th element of the inverse of the N x N correlation 
matrix C that characterizes the multivariate Gaussian. We refer the interested 
reader to ref. 



for a derivation of (2.121 and (2.131 



2.2. Statistical Symmetries 

From the observation of turbulent flows it is immediately clear that individual 
realizations of turbulent fields break spatio-temporal symmetries. Depending 
on the flow under consideration, however, symmetries may be recovered in a 
statistical sense [37], i.e. due to the averaging procedure. The rigorous study 
of statistical symmetries has been introduced to hydrodynamical turbulence 
by Robertson [3S] , Chandrasekhar |39l HP] and Batchelor [IT] and since then 
has been an indispensable mathematical tool for the statistical investigation of 
turbulent flows. The Howarth-Karman relation [32] may serve as a prominent 
example here demonstrating that statistical symmetries may not only help 
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to simplify the statistical equations, but also may cast them into a tractable 
mathematical form allowing for new insights. Of course this is also the case for 
the application of statistical symmetries to kinetic equations, and in fact, some 
of the results presented in the following depend crucially on the exploitation of 
these symmetries. 

Among the simplest of these symmetries are statistical stationarity, homogene- 
ity and isotropy, and each of these symmetries can be defined as the invariancc 
of a given statistical quantity with respect to a symmetry transformation. In 
case of stationarity and homogeneity, the corresponding transformations are 
temporal and spatial shifts, respectively. For the TV-point vorticity PDF of 
two-dimensional turbulence this leads to the following relations 

/({«,,*,}, t) = f({u h x x + r}, t) (2.14a) 
/({«!, xi}, t) = f({0J hXl }, t + r), (2.14b) 

where r and r are arbitrary temporal and spatial shifts. Isotropy is defined as 
invariance with respect to rotations and (optionally) reflections. For the iV-point 
velocity PDF this leads to 

f({u h xi}, t) = f({Ru h Rx t }, t) , (2.15) 

where R is an arbitrary rotation matrix. For the vorticity PDF one has to take 
into account that the vorticity is a pseudo-scalar, i. e. it changes sign when R 
contains a reflection. 

To demonstrate the importance of these symmetries let us consider the 
single-point velocity PDF f(u,x,t) in homogeneous, isotropic and stationary 
three-dimensional turbulence. As the PDF depends only on one spatial and one 
temporal point, the invariance under the corresponding shifts implies indepen- 
dence of these variables, i.e. f(u,x,t) = f(u). The invariance under rotations, 
i. e. f(Ru) = /(it), implies that the PDF can only depend on the magnitude of 
the velocity vector instead of the three vector components. Counting all of these 
simplifications together, we have reduced the dependency of the single-point 
PDF from seven (one temporal, three spatial and three sample-space) variables 
to one, the velocity magnitude. 

The fact that a PDF does not depend on certain variables (due to statistical 
symmetries) allows one to formulate a whole sequence of relations which involve 
conditional expectations. Let us exemplify this for the PDF of a field u>(x,t) 
with stationary statistics. The time derivative as well as the higher-order time 
derivatives have to vanish identically. This leads to 



d . d 

0=-/(u,,*)=-- 



m u(x,t) 



did, . 



S(lj — lu(x, i)) 
u\f(u,x). (2.16) 



Here we already encounter the closure problem, from eq. (2.161 it is clear that 



we need the joined PDF of Mru: and uj to determine the temporal variation of 
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f(uj,x,t). The resulting equation for £^f(cj,x) 
Pope and Ching [26]: 



has been highlighted by 







did 2 , . 



uj )f(u,x) + 



du; 2 



uj(x, t) 



f(u,x) (2.17) 



Obviously this treatment can be generalized to all continuous symmetries. For 
translationally invariant statistics the simplest relation reads 







d_ 

duj 



[Au(x,t)\u)f(u,t) 



duj- 



■([Vu(x,t)} 2 \u)f(u;,t) 



(2.18) 



This type of relation will play an important role in the treatment of dissipative 



terms in the evolution equations, see section 4.2 



Up to now we have applied the symmetries to PDFs, which are scalar-valued 
functions. However, the kinetic equations of the LMN hierarchy also contain 
vector- and tensor-valued functions, e.g. (Au\u). Their representation can also 
be simplified with the help of statistical symmetries, which leads to the general 
concept known as the invariant theory of tensors. Generally speaking, exploiting 
statistical symmetries, vectors and tensors can be expressed as a superposition 
of a generating set of covariant vectors or tensors. 

Here we would like to exemplify this concept for isotropic vector- and tensor- 
valued functions depending on one vector variable. Such functions appear in 
the context of the LMN hierarchy as conditional averages of vectors and tensors 
conditioned on a vector, e.g. (Au\u) and (Vu(Vu) T |u). Due to the fact that 
in two dimensions fewer reductions are possible we restrict the further discussion 
to three dimensions. 

Let (a | it) and (Bij\u) be a conditionally averaged vector and symmetric 
tensor in three dimensions. Statistical isotropy demands that they obey the 
transformation rules 



(a\Ru) = (Ra\u) 
[B\Ru) = (RBR T |u) , 



(2.19a) 
(2.19b) 



which can be derived from the transformation rule of a PDF (2.15). Due to 



these symmetries this conditional averages must have the representations 



'B i:j \u) 
X(u) 



(u ■ a\u)u 

n{u)5ij + [A(u) - n{u)]uiUj 
(u ■ Bulu) 



^Tr(B)|u) 



-(« ■ Bu\u) , 



(2.20a) 
(2.20b) 
(2.20c) 

(2.20d) 



where u = u/u is the unit vector in direction of u. These results can be rigorously 
derived, where the core idea is to use rotations around the axis u. 

These representations have several advantages. On the one hand they make 
it possible to simplify the equations containing this conditional averages. On the 
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other hand the symmetry reduction renders the estimation of the conditional 
averages from DNS feasible as only scalars conditioned on scalars have to be 
estimated from the data, in contrast to vectors or tensors conditioned on vectors. 

Finally, we would like to mention here that statistical symmetries can be 
used to formulate further constraints on the conditional averages like (u ■ a\u) 
and X(u). For example, A(u) has to be an even function in u, i.e. its Taylor 
series may contain only even powers of u. For further details we refer to [43] , 

2.3. LMN Hierarchy for the Vorticity in Two-Dimensional Turbulence 

We now proceed to a detailed presentation of the hierarchy of evolution 
equations for the PDFs. To keep things simple, we focus on the hierarchy for 
the vorticity field in two-dimensional turbulence. The corresponding, analogous 
hierarchies for the three-dimensional case will be summarized later. 

We start from the evolution equation for vorticity ui = Lue z , where the 
vorticity component u>(x,t) depends on the horizontal coordinates x — (xi,^). 
The basic fluid dynamic equation for the vorticity reads 

d 

t) + u(x, t) ■ V x u(x, t) = L(-A x )uj(x, t) + F{x, t) . (2.21) 

We have included a forcing term, F(x,t) as well as a general viscous term, 
L(—A)u>(x,t), which may contain large-scale friction, usual dissipation and 
hyperviscous terms. A suitable choice of F and L allows one to generate a direct, 
an inverse |36J as well as dual cascades (see e.g. [44]). The velocity field can 
be determined using Biot-Savart's law, which relates the vorticity field to the 
velocity field: 

u(x,t) = [dx'uj(x',t)e z x — X - , = [dx'U(x - x f )u(x',t) (2.22) 
J 2tt\x-x'\ 2 J 

The velocity field U(x — x') may be regarded as the field of a single point vortex 
located at x' . It is purely azimuthal and decays like l/|cc — x'\. 

We now want to proceed to a statistical description of the vorticity. To derive 
the kinetic equation for the single-point statistics of two-dimensional turbulence, 
we can use the concepts introduced in section [2] in a straightforward manner. To 
this end we consider the fine-grained PDF 

A(wi,xi,t) = 5(wi-w(xi,t)) (2.23) 
and take the time derivative, which leads to 

d 



^/i(wi,xi,t) = 



OLOl 



gt u(x u t). (2.24) 



Analogously we compute 

ll(xi,t) ■ Vas,. A ("1 



u(x u t)- V a!l a;(xi,i). (2.25) 
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Adding these results we obtain 
d 



dt 



fi(ui,xi,t) + u(x 1 ,t) ■V Xl fi(uJi,X 1 ,t) 



d_ 
m 



/i(wi,xi,t) + V Xl ■ [tt(aci,t)/x(wi,xi,t)] 



d 
duji 
d 

d 



fi(uii,xi,t) 

fl{u)i,Xi,t) 



u>(xi,t) + u(xi,t) ■ V Xl u>(xi,t) 



[L{-A ai )u(x lt t) + F{x u t)} 



[L{-A Xl )u{ Xl ,t) + F(a;i,t)]/i(wi, xi, t) . 



(2.26) 
(2.27) 
(2.28) 

(2.29) 
(2.30) 



For eq. (2.271 we have used solenoidality of the velocity field. From eq. (2.28) 
to eq. (2.291 we have replaced the left-hand side of the vorticity equation (2.21 ) 
with its right-hand side. The step from eq. (2.29) to eq. (2.301 finally makes 
use of the fact that the terms of the vorticity equation are independent of the 
sample-space variable, which concludes the derivation of the single-point equation 
for the fine-grained PDF. Ensemble averaging then leads to the desired PDF 
equation. 



dt 



fl(u]l,Xi,t) + Va,! • (u(x 1 ,t)f 1 (uj 1 ,X 1 ,t)) = 



(2.31) 



d 
duji 



{[L(-A ai )u(xi,t) + Fixt^MuuXi,*)) , (2.32) 



which confronts us with the closure problem in terms of the joint averages of the 
fine-grained PDF and velocity, the dissipative term and the forcing term. To 
proceed, two options are available. First, the joint averages can be expressed as 
conditional averages according to 

(u(x 1 ,t)f 1 (u 1 ,x 1 ,t)) = (u(x 1 ,t)\uj 1 ,x 1 ,t)f 1 (uj 1 , Xl ,t) (2.33) 

([L(-A Xl )uj(x 1 ,t)]f 1 (uj 1 ,x 1 ,t)) = (L(-A Xl )uj(x 1 ,t)\uj 1 ,x 1 ,t)f 1 (uj 1 ,x 1 ,t) 

(2.34) 

= (F(xi, t)|u*,xi, a:i,t). (2.35) 

These terms are in principle accessible by DNS data, which will be discussed 
in detail below. The second option is to express the unclosed velocity and 
dissipative terms by the two-point statistics. For the velocity term we obtain by 
insertion of the identity J dw2 5(ui2 — oj(x2,t)) = 1 



(u(xi,t)f(ui,xi,t)) = J dx 2 J du 2 U(xi - x 2 ) uj 2 f(ui,xi,u} 2 ,x 2 ,t) . (2.36) 

At this stage the coupling to the two-point PDF arises, in much the same way 
as for the Friedmann-Keller hierarchy for the correlation functions. A similar 
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treatment is then performed for the dissipative term, which leads to 
([if-A^Mxi, *)]/(«!, Si,t)) 

= jdx 2 Jduj 2 5(x 2 - xi) u 2 L(-Ax 2 )f{vi,Xi,L)2,X2,t) . 

Up to now, we have not further specified the forcing term. For a deterministic 
forcing term, the unclosed term simply involves the joint statistics of the external 
force field and the vorticity field. Especially for two-dimensional turbulence a 
white- in-time forcing term with Gaussian statistics is often considered, which 
can be treated further. Such a forcing term is fully specified by the correlation 
function 

(F(x 1 ,t)F(x 2 ,t / )) = Q( Xl - x 2 )6(t - t') . (2.37) 

The transition from the fine-grained PDF to the full PDF then involves an 
averaging procedure with respect to the white noise force field F(x,t). This 
yields a representation of the forcing terms as 

1 2 d 
(F(x 1 ,t)f(uj 1 ,x 1 ,u; 2 ,x 2 ,t)) = --^2Q(xi - x j ) — f(u; 1 ,x 1 ,u: 2 ,x 2 ,t) . 

(2.38) 

The derivation of this result is similar to the derivation of the Fokker-Planck 
equation for a set of Langevin equations |451 146] . The random force method 
applied to turbulence has also been discussed by Novikov (47]. 

The whole derivation is generalized to the TV-point statistics in a straight- 
forward manner by starting from the TV-point fine-grained PDF. Expressing the 
unclosed terms in terms of the (N + l)-statistics then results in 




where again the coupling to the (N + l)-point PDFs arises due to the nonlocal 
Biot-Savart law relating velocity and vorticity field, as well as the dissipative 
term. The alternative formulation is obtained by replacing the terms containing 
the (7V + l)-point probability functions by the corresponding conditional averages. 
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The set of kinetic equations then takes the form 



a 

dt 



f({u>i,xi},t) + v ^ • (M(^i,i)l{w/,^;})/({w;,a; ; },i) 



3=1 

N 



= - E ^M^,*)|{^}>/({^ *) . (2-40) 

3=1 3 

where we have introduced the abbreviation for the right-hand side of the vorticity 



equation (2.21 1 



n(x J ,t)=L(-A mj )u(x J ,t) + F(x j ,t). (2.41) 

In this formulation, the chain of PDF equations is retained up to the order N 
and formally truncated on this level by the introduction of the up to this point 
unspecified conditional averages. The case of N — 2 will be discussed with the 
input from DNS data in section [3J 

2.4- The Method of Characteristics 

The kinetic equations formulated in terms of conditional expectations formally 
are first-order partial differential equations. Hence, they can be solved using the 
method of characteristics [48] . By application of this method, one can identify 
trajectories in the phase space spanned by the variables of the partial differential 
equation, the so-called characteristic curves or simply characteristics. Along 
these characteristics, the original partial differential equation transforms into an 
ordinary differential equation. 

Let us exemplify the procedure for the simplest case of the kinetic equations for 



two-dimensional flows, i.e. the ./V-point kinetic equation (2.40). The associated 
characteristics are solutions of the equations 

-f-t(s) = 1 (2.42a) 
as 

= (M(^j^)l{ w ^ a; i}>{^,x i }={f2 i ( s ),x i ( s )} (2.42b) 

-^%(s) = (M(a5i^)|{wi,Xi}){ tl ; I , a!l }={n J (»),x,( a )} > (2.42c) 

where s is the parametrization of the trajectories through phase space, and the 
conditional averages have to be evaluated at the current positions along the 
characteristics. The solutions of these ordinary differential equations are the 
characteristics Xj(s, {to®, x°}) and fXj(s, x°}) that are depending on the 
initial conditions {ojf,x®} at s = 0: 

X J (0,{u J °,x°})=x° (2.43a) 
Q j (0,{ulx°})= C j° (2.43b) 
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Along these characteristics, the kinetic equation (2.40) transforms into an ordi- 
nary differential equation: 

d N 

^/OO = -E[ V -, ' (tt(a?i,t)|{wi,x,}> 

3=1 

"Xxj-.tJKw,,*,})! /(s) (2.44) 



{ldi,:E(} = {fi;(s),X;(s)} 



This means that the deformation of the PDF along the characteristics is governed 
by the bracketed term on the right-hand side, which is the phase space divergence 
of the vector field formed by the conditional averages appearing in eq. ( 2.42| . 



If the PDF at the initial point {uf , x®} of the characteristic curve, /o({ w / U : xf}), 
is known, the PDF at "later" points Xj = Xj(s, {u)f, x?}), uij = f2j(s, {w; , x®}) 
(i.e. for a certain value of s > 0) is determined by 

f{{u l ,x l }) = J/ (K,x°}). (2.45) 

Here, J denotes the Jacobian determinant of the inverse mapping (X J , QJ 1 ) 
that matches the phase space points Xj and ljj back to the initial conditions of 
the characteristic that passes through them, i. e. 

Xr 1 ({x l ,ui}) = x° j (2.46a) 
nj 1 ({xi,u l }) = u)?. (2.46b) 

The method of characteristics exhibits some interesting analogies to the 
Lagrangian description of turbulence. To this end one can note that evolu- 



tion equations of the characteristic curves, eqs. (2.42), describe the motion of 
Lagrangian fluid particles that are moving within the conditionally averaged 
velocity field and change their vorticity according to the conditionally averaged 
dissipation and forcing fields. In this sense the characteristic equations describe 
the average dynamics of a class of fluid particles defined by the configuration 
of the N fluid particles under consideration. This view is closely related to the 
concept of quasi-particles. 

Furthermore, in the continuum limit N — > oo or (Xj , CljJ — > (X(y, t), Q(y, t)j , 
where the characteristics are identified by their continuous starting positions y 
that cover the whole space, the characteristic equations tend to the vorticity 
equation formulated in the Lagrangian framework, 

±X(y, t) = J Ay' n(y', t) U(X(y, t) - X(y', tj) (2.47a) 

t) - [£(-A)w(x, t) + F(x, t)] x=x{y t) , (2.47b) 

with U again being the Biot-Savart kernel. The formulation, modeling, and 
investigation of the characteristic equations, therefore, allows one to make contact 
with models for the motion of few Lagrangian tracers, especially tetrad models 
|491 1ST)] and models for the velocity gradient tensor |511 152"] . We refer to the 
reviews of Pumir and Naso, and Li in this volume. 
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2.5. From the Eulerian to the Lagrangian Frame 

The discussion of the preceding section links our discussion of the LMN 
hierarchy to the Lagrangian description of turbulence. An extended discussion of 
the Lagrangian approach would be beyond the scope of this review paper but we 
want to make at least some comments on this issue. The choice of the frame of 
reference depends mainly on the physical problem under investigation. Whereas 
the Eulerian picture seems for example preferable for questions connected to 
the energy flux in scale, the Lagrangian picture is the natural way to describe 
transport and mixing processes. To this end one considers the Lagrangian 
path. X(y,t). of a Lagrangian particle starting at initial point y at time t = 0. 
Furthermore, one introduces the Lagrangian velocity field or vorticity field 
according to u^(y,t) = u(X(y,t),t), u}^(y,t) = u>(X(y,t),t). The Lagrangian 
single-point PDF is then defined as 

f L (x, u, y, t) = (S(x - X(y, t))S(u - u L (y, t))) (2.48) 

with a straightforward extension to N particles. As for the Eulerian case one 
can now proceed to construct a LMN hierarchy for the Lagrangian quantities of 
interest. For example the chain of equations for the case of the velocity PDF 
has been formulated in |53l 153] , whereas the case of the vorticity PDF has been 
investigated in [55]. 

Another important topic is connected to the general relation between the 
Lagrangian and the Eulerian picture. In some situations it is helpful to under- 
stand the statistics of an observable in one frame of reference given the statistics 
of this observable in the other frame of reference. One famous example is the 
translation of the two-point two-time velocity autocorrelation function from 
the Eulerian into the Lagrangian picture (see e. g. j5H] and references therein). 
In recent years especially the statistics of Eulerian and Lagrangian velocity 
increments have gained much interest. Experiments and numerical simulations 
showed that the deviation from Gaussianity of these PDFs is in the Lagrangian 
frame even stronger than in the Eulerian frame. Attempts to relate increment 
PDFs and structure functions have been made in Refs. [57] [55] [55] mainly based 
on dimensional arguments. 

In |60| the formalism of fine-grained distributions introduced in section [2] was 
used to derive an exact relation of these quantities which includes the transition 
relations used in [57] [55[ [59] as a special case. The main result of [50] is that in 
general one has to add information in form of transition PDFs to translate the 
PDFs of an observable from one frame of reference to the other. The resulting 
relation that translates the increment PDFs involves two transition PDFs which 
add the missing information. One takes into account the dispersion of the 
Lagrangian tracers and the other adds information on the time development of 
the flow field at the starting point of the tracer. With the help of this relation it 
was for example possible to understand how the nearly Gaussian PDFs of the 
Eulerian velocity increments in the inverse cascade of two-dimensional turbulence 
translate into non-Gaussian PDFs in the Lagrangian frame. An application to 
three-dimensional hydro- and magneto-hydrodynamic turbulence can be found 
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Figure 1: Visualization of the vorticity field in two-dimensional turbulence in the inverse cascade 
regime. The vorticity field is taken from simulation A in |36| . Videos of this visualization can 
be found in 1621 . 



in |61J. In principle the methods to derive the results of [60 can be used in 
general to find exact relations between Eulerian and Lagrangian PDFs for various 
observables. 



3. Inverse Cascade in Two-dimensional Turbulence 

In this section we review results of an analysis of two-dimensional turbulence 
based on the kinetic equation approach and the results derived in section [2. 3 1 
The performed analysis along the lines of the LMN hierarchy is based on a DNS 
of two-dimensional turbulence in the inverse cascade regime. The DNS data 
comes from a pseudospectral simulation of the two-dimensional Navier-Stokes 
equation in vorticity formulation (cf. fig.[lj. Forcing at small scales generates 
an inverse cascade, where energy is transported to large-scale motions, which 
are destroyed by large-scale dissipation in order to generate stationary statistics. 
These simulations show a constant energy flux within the inertial range. Details 
on the setup of the simulations can be found in |36j . In the following we present 
the conditional averages appearing at the second level of the LMN hierarchy 
(obtained from DNS) and compare them with a Gaussian approximation. 



3.1. Two-Point Kinetic Equation 

The properties of the vorticity field in the inverse cascade has been analyzed 
on the basis of the two-point statistics. Considering eq. (2.401 for the case of 
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Figure 2: Conditional vorticity field for two fixed vorticities in two-dimensional turbulence in 
the inverse cascade regime. The vorticity field is localized near the two base points and shows 
similarities to the vorticity field of two vortex patches. The induced velocity field is marked 
with vectors and streamlines. Left: Measu red d irectly from DNS. Right: Re-constructed 
following the Gaussian approximation in eq. 



N = 2 and additionally assuming stationary statistics, we obtain 

[V a . 1 -{u(x 1 ,t)\L) 1 ,X 1 ,LJ2,X2) + V X2 • (u(x2,t)\ui,X 1 ,LJ2,X2)]f(ui,X 1 ,U2,X 2 ) 

- d 

{(X(xi,t) \iVi,Xi,Ul2,X 2 ) + 

-ouji 

d 1 

- (^(x 2 ,t)\uj 1 ,X 1 ,UJ2,X2) f(uJ 1 ,X 1 ,UJ 2 ,X 2 ). (3.1) 

0U)2 J 

As has been shown in [35], the unclosed terms are specified by the conditional 
vorticity field {uj{x, t)|wi, X\, Ld2, x 2 ), which especially allows us to evaluate the 
conditional velocity field (u(x,t)\uji,xi,uj2,x 2 ). In |36) these two terms have 
been determined from DNS data. The result is shown in fig. [2] The fields 
reveal close similarities with fields generated by a vorticity field which is strongly 
localized at the base points xi, x 2 . Therefore, the characteristic equations 
associated to the kinetic equations for few points are closely related to the 
dynamics of localized vortex patches, which yields a direct connection between 
the statistical description of two-dimensional turbulence and vorticity dynamics. 
In the simplest case these patches can be approximated by point vortices. Hence, 
the knowledge about the dynamics of point vortex systems (we refer the reader 
to the monograph of Newton [63 and the review articles by Aref [53]) turns out 
to be of great importance. Furthermore, the conditional dissipation and forcing 
field turns out to be a linear function in wi, uj 2 , see |36J. As a consequence 
one may visualize the two-point statistics on the basis of vortex models for two 
localized vorticity patches. 

3.2. Gaussian Approximation to the Conditional Velocity Field 

The simplicity of the above observation has motivated further analytical 
investigation in [36]. A first approximation to the conditional velocity field 
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can be obtained on the basis of a Gaussian approximation of the conditional 



PDF. The general expression (2.121 evaluated for Gaussian statistics yields the 
following expression for the conditional vorticity field: 

(uj(x,t)\uj 1 ,Xi,UJ2,X2) 

= C(x - Xl ) [C^ui + C^w 2 ) + C{x - x 2 ) [C^V + C^-wz] (3.2) 
Here, C{x — x') = (lj(x, t)u>(x', t)} is the two-point correlation function of the 



vorticity field (cf. eq. (2.131 for the definitions of the further quantities). This 
revealing result shows that the conditional vorticity field in fact is made of vortex 
patches with vortex profiles given by the two-point correlation. 

The conditional velocity field can be determined from an application of 



Biot-Savart's law (cf. eq. (2.22 1). The evaluation of the corresponding formula 



leads to the introduction of a screened velocity field: 

U s (x-Xi)= fdx'U(x-x')C(x' -Xi) (3.3) 

The subscript S shall indicate that the velocity field Us may be regarded as 
the field of a screened vortex similar to the notion of Landau quasi-particles. 
For a point vortex, the velocity field is obtained by replacing C(x — x') with a 
(5-function. For widely separated base points the conditional velocity field turns 
out to be a superposition of the fields of two vortices localized at X\, x 2 : 

(u(x,t)\LJx,X 1 ,U2,X 2 ) = -^TF7^\Us{x ~ X\)u\ + U S (x - X 2 )bJ 2 ] (3.4) 

However, the velocity profile is different from the one of point vortices, due to the 
fact that the turbulent surrounding has been included in terms of the conditional 
Gaussian approximation. This fact allows one to visualize the characteristic 
equations as describing the dynamics of quasi-vortices in the sense of Landau 
quasi-particles. Fig. [2] compares the conditional velocity fields determined from 
DNS with the Gaussian approximation. 

3.3. Failure of the Gaussian Approximation 

Although the Gaussian approximation yields a reasonable representation 
of the conditional vorticity and velocity field it does not reproduce the energy 
transport of the inverse cascade. The conditional velocity fields 

(u{x u t)\u x ,xi,u) 2 ,x 2 ) = cJ^ U s( x i ~ ^2)^2 
(u(x 2 ,t)\uJi,xi,u 2 ,x 2 ) = cJojUsix? - (3.5) 

only have a component in the direction e 2 x (xx — x 2 ), which has important 
implications for the two-point statistics. This can be seen by noticing that for 
homogeneous flows the left-hand side of eq. ( |3.1[ ) can be rewritten as 

[V Xl • (u(x 1 ,t)\uj 1 ,x ll uj 2 ,x 2 ) + V X2 ■ (u(x 2 ,t)\uj 1 ,x 1 ,uj 2 ,x 2 )}f(uj 1 ,x 1 ,uj 2 ,x 2 ) 
= [V r • (u(x 2 ,t) - u(x 1 ,t)\uj 1 ,x 1 ,uj 2 ,x 2 )f(uj 1 ,x 1 ,uj 2 ,x 2 )} (3.6) 
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with r = x 2 — x\. As a result, this convective term vanishes for the Gaussian 



approximation. Calculating moments from eq. (3.1 1 shows that it is precisely 
this term which is related to the enstrophy and energy fluxes across scales, which 
then also implies that these fluxes vanish in a Gaussian approximation. For a 
more detailed discussion of the Gaussian approximation we refer the reader to 

A particularly intuitive interpretation of this statistical result is achieved 
with the help of the method of characteristics, which evolve according to 

Xl = W) s{xi ~ Xs) 

° Jl TT I \ 

X2 = 6(0) s ^ 2 ~ Xl > 

LOl = (tl(xi,t)\u>i,Xi,U)2,X 2 ) 

6j 2 = dj,(x 2 ,t)\u) 1 ,x 1 ,u)2,x 2 ) ■ (3.7) 

This set equations shows that statistically an azimuthal motion with time- 
dependent vorticity amplitudes is induced, a change of distance however is 
absent . 

What apparently is missing in the Gaussian approximation is a relative 
motion of the two quasi- vortices, i.e. a motion in the direction of r. Therefore, 
a straightforward generalization of the characteristic equations takes the form 

W2 TTt \ , X l~ X 2 rjn I \ 

= TT77^ u \ x i - X V + ] ; H (\ x i - X2\,ui,U2) 

6(0) \xi-x 2 \ 

X 2 = 7^Y C/ ( a; 2 - Xi) + —^77(1^2 - Xl\,U) 2 ,Ul) , (3.8) 

6(0) \x 2 -xi\ 

with a function H that depends on distance only. This relative motion is the 
signature of the inverse cascade in the characteristic equations associated to the 
kinetic two-point equation: The surrounding fluid turbulence leads to an effective 



attraction or repulsion of the quasi- vortices in the characteristic equations (3.8). 

A successful theory of the inverse cascade has to provide an answer to the 
question on the dynamical origin of the relative motion of two quasi- vortices. 
Recently, it has been suggested that the relative motion is induced by an elliptical 
deformation of the vorticity profile of the vortices, due to the mechanism of 
vortex thinning, i. e. an elongation and thinning of the quasi-vortices due to its 
mutual interaction. Modeling these effects on the basis of a generalization of 
Onsager's point vortex models by introducing elliptical point vortices which can 
respond to the large-scale field generated by other vortices has been shown to 
be successful. The corresponding dynamics exhibits the existence of an inverse 
cascade by the formation of clusters of vortices with the same circulations |65j . 
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4. Velocity Statistics in Three-Dimensional Turbulence 



4--1- LMN Hierarchy for the Velocity 

The hierarchy of evolution equations for the iV-point probability distributions 
of the velocity field 



f({u l ,x l },t) = {l[S(u l -u(x l ,t)) ) (4.1) 



\Z=1 / 

has been formulated independently by Monin and Lundgren [TOlUT]. The starting 
point for the derivation is the Navier-Stokes equation for an infinitely extended 
incompressible fluid 

^+u{x,t)-v\u(x,t) 
= -Vp + vAu(x, t) + F(x, t) 

= - J dx' K(x - x')V x , ■ [u(x', t) ■ V x ,u{x\ t)] + vAu(x, t) + F(x, t) , 

(4.2) 

where the pressure term has been explicitly represented as a functional of the 
velocity field using the kernel 

K{x-x') = ^^^- s . (4.3) 



■in 



x — x' 



The field F(x, t) is an external force which drives the fluid motion. For the sake 
of simplicity we neglect the forcing in this section. The evolution equation for 
the iV-point velocity PDF can be obtained with the method of fine-grained PDFs 
as discussed before, resulting in 

d N \ 

N 

^V Mj ■ du'dx'K(x 3 -x')(u' ■V x , ff(u',x',{u h x l },t) 

3 = 1 J 

N 

-v22 Vu i' du'dx'S(x j -x')u'A x >f(u',x',{u u xi},t). (4.4) 

3=1 

A readable discussion of the steps involved in the derivation can be found in 
ref. [TD]. This equation couples to the (N + l)-point PDF through the two terms 
on the right-hand side, which represent the contributions due to the pressure 
gradient and the Laplacian of the velocity field. Therefore what we just obtained 
is an infinite hierarchy of coupled evolution equations, the LMN hierarchy for 
the velocity. 
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As a side note, we want to mention that the Friedmann-Keller equation 
for the correlation functions (ui(x\,t) . . .Uj(xj\r,t)) can be obtained from the 
LMN hierarchy by representing the correlation functions as moments of the 



corresponding PDFs and using the evolution equation (4.4). 

As discussed before, it is possible to truncate the hierarchy at any level 
through the introduction of conditional averages for the pressure gradient and 
the dissipation. Doing this at the iV-th level we obtain 



dt~ 



N 



• V Kj .)f{{ui,xi},t) = 
4=1 ' 

N 

4=1 
JV 

-^V Uj • {v/^ Xj u(x h t)\{u u xi})f({u u xi},t) 
4=1 

JV 

-^V Uj ■(F(x j ,t)\{u l ,x l })f({u l> x l },t), (4.5) 
4=1 

where we have now included the forcing. This evolution equation is formally 
closed but contains the unknown conditional averages of the pressure gradient 
and the dissipation. These terms can be either modeled or estimated from 
experiments or simulations. 

4-. 2. Single-Point Velocity Statistics 

We now come to a closer investigation of the single-point velocity PDF in 
homogeneous isotropic turbulence. Although this is maybe one of the simplest 
statistical quantities to consider, there has been an ongoing discussion in the 
past whether this quantity exhibits deviations from Gaussianity. Gaussianity is 
often assumed based on an argumentation exploiting the central limit theorem. 
Numerical and experimental investigations, however, give evidence for slight 
deviations from Gaussianity, especially indicating the occurrence of sub-Gaussian 
tails [66, 67, 68J. From the theory side, varying opinions can be found in the 
literature. Falkovich and Lebedev have used the instanton formalism to argue 
in favor of sub-Gaussian tails depending on the external forcing [69J . Gaussian 
PDFs have been found for decaying turbulence by Ulinich and Lyubimov |15] 
and later by Hosokawa |2Q] or in the case of the cross-independence hypothesis 
by Tatsumi and coworkers |21| . 

By combining the kinetic equation of the single-point velocity with DNS 
results, we will argue here for the existence of deviations from Gaussianity caused 
by statistical correlations, e. g. of the dissipation field and the velocity field. For 
more details we kindly refer the reader to refs. |4"5 t ITU j ITT]. 

From now on we will drop the argument (x, t) for all realizations. We remind 
the reader, however, that the first argument of a conditional average is understood 
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as a realization whereas the condition is specified by the corresponding sample- 
space variable. 

Introducing conditional averages on the level of the first equation of the 
hierarchy we arrive at 



d_ 

at 



f(u,x,t) + u-'Vf(u,x,t) = -V«- (-Vp + vAu + F\u) f(u,x,t) . (4.6) 



The left-hand side is closed, however, the right-hand side contains unclosed terms 
involving the joint statistics of the pressure gradient, viscous diffusion and the 
external forcing with the velocity. 

To make a connection to DNS results, this equation has to be simplified with 
the help of statistical symmetries. As exemplified above, the PDF is independent 
of the spatial coordinate for homogeneous flows. As a result, the advective term 
on the left-hand side vanishes. Homogeneity further lets us derive the relation 



BP 
dx 



2 







_d_ 
dui 



d 2 Ui 



u )/(«,*) + 



d d I dui duj 



dus, dv,j \ dxi. dx 



u)f(u,t) 



which is a multi-dimensional generalization of the relations introduced and 
investigated by Ching and Pope [26, 2S][27|. Inserting this result into the kinetic 
equation (4.6 1 leads to 



J)_ 



O.i 



u)f(u,t) 



duj du 



i L L 

dx k dx k 



«^ /(«,*)■ 
' . . (4 ' 8) 

By this we have replaced the conditional diffusive term with the conditional 
(pseudo-)dissipation tensor 



Dij{u) 



dui du-j 
i ± 

dxh dxi 



u 



(4.9) 



For isotropic turbulence, all of the appearing quantities have to be invariant 
with respect to rotations. As a consequence, the velocity PDF depends on the 
velocity magnitude only and is determined by the PDF of the velocity magnitude 
by the relation 

f(u,t)=4nu 2 f(u,t). (4.10) 



The conditional vectors appearing in (4.6 1 take the form 



(— Vp|it) = n(u)u 



(F\u) 



A(u)u 
$(m)m 



n(u) = (—u ■ Vp\u) 

A(u) — (vu ■ Am|m) 
<i>(u) = (u ■ F\u) . 



(4.11a) 
(4.11b) 
(4.11c) 



Note that the conditional averages II, A and $ are scalar functions of a scalar 
argument, which makes an estimation from DNS data feasible. The conditional 
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dissipation tensor can be simplified in a similar manner and can be written down 
in terms of its eigenvalues according to 

Aj(u) = n(u)S lJ + [A(u) - (4.12a) 

fi(u) = i(e + ^ 2 |w) - i(^(A T u) 2 |u) (4.12b) 
A(u) = (v{A T u) 2 \u) . (4.12c) 

Here, the conditionally averaged kinetic energy dissipation, the squared vorticity 
and the transpose of the velocity gradient tensor Aij = enter, clarifying that 
the single-point velocity statistics depends on the correlation of the dissipation 
field and the velocity field. Taking all of these results together, eqs. (4.6) and 
(4.8 1 can be recast in the forms 

at /(m ' t] = ~ du ( n(u ' t] + A(u ' t] + ^ t] (4 ' 13a) 

~f(u, t) = ~ (n(«, t) + *(«, t) - 2/ ^ ) ) /(«, t) - ~A( W) *)/(«, t) • 

ot ou\ u J ou z 

(4.13b) 

Two interesting conclusions can be drawn for stationary flows. First, eq. ( 4.13a| 
implies 

U(u) +A(u) + $(m) = 0, (4.14) 

i. e. balance of the pressure term, the diffusive term and the forcing term for 
each fixed value of velocity. Even more importantly, a stationary solution can 
be derived for (4.13b I taking the form 

f M = TT~\ exp / du ' 4 - 15 

where Af is a normalization constant depending on u . This solution is an exact, 
yet unclosed relation. Without further input of numerical data, the conclusion 
to be drawn from these theoretical results is that the shape of the PDF is 
determined by the conditional averages which encode the correlations of the 
various turbulent fields appearing in the Navier-Stokes equation with the velocity 
field. 

Before turning to the DNS data, we would like to demonstrate how a "naive", 
yet rational closure can be established by simple physical arguments [70]. To 
this end we make use of the fact that the conditional averages can be reduced to 
ordinary averages by integration, which lets us introduce a number of integral 
constraints. For example, for homogeneous isotropic flow, the pressure term has 
to vanish, implying 

/•oo poo 

= (u-Vp)= du(u- Vp|u)/(u) = - / duull(u)f(u). (4.16) 
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Figure 3: Visualization of the velocity (left) and vorticity field (right) in three-dimensional 
homogeneous isotropic turbulence. The figures show the color-coded magnitude of the full 
vector, where blue represents low and red high values. Additional visualizations to be found at 

El. 



The dissipative terms have to yield the mean kinetic energy dissipation 

poo 

-(e) = (vu- Att) = / duuA{u)f(u) (4.17a) 







(e) = (Tr(D)) = / du [X(u) + 2(i(u)]f(u) . (4.17b) 



n 



As for stationary flows the forcing terms balances the mean energy dissipation, 
we obtain 

/■oo 

(e) = (u-F)= duu$(u)f(u). (4.18) 
Jo 

Based on these constraints, a lowest-order approximation of the conditional 
averages yields 

n (u) = (4.19a) 



$o(u) = (4.19b) 
so 1 



Ao(u) = -^«, (4.19c) 

where a = W (u 2 /3). In this approximation, pressure effects are neglected 
completely and the forcing and diffusive term are linear functions of the sample- 
space velocity. For the conditional dissipation tensor we arrive at 

Ao(w) = no(u) = -y- , (4.20) 
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i. e., the two eigenvalues are identical and can be expressed in terms of the average 
rate of kinetic energy dissipation. This last approximation can also be understood 
in terms of a simple decoupling argument: Assuming a scale-separation between 
the small-scale dissipation field and the comparably large-scale velocity field 
(see fig. [3] for a visualization of the velocity and the vorticity field) implies that 
the conditional averages reduce to ordinary averages. Evaluating the stationary 
solution ( |4.15 1 with these approximation yields a Maxwellian distribution of the 



velocity magnitude 

l~2u 2 ( u 2 \ 

^=v^ e H"^J' (4 - 21) 

which corresponds to a Gaussian distribution of the velocity vector, 

1 / u 2 \ 

f^ = j2^W^{-^)- (4 ' 22) 

This interesting result shows that a Gaussian single-point velocity PDF can be 
obtained as a lowest-order approximation. The DNS results in the following, 
however, will show that pronounced deviations from this approximation exist 
and in combination lead to deviations from Gaussianity. 

The DNS results have been obtained with a standard pseudospectral simula- 
tion of the vorticity equation with periodic boundary conditions. The forcing is 
acting on large scales and rescales the amplitudes of the Fourier coefficients in a 
wavenumber band such that the total energy of the flow is conserved while letting 
their phases evolve freely. The simulation exhibits a Taylor-based Reynolds 
number of R\ m 112 at a resolution of 512 3 grid points. To obtain good ensemble 
statistics, the simulation was run for more than 150 large-eddy turnover times 
and two hundred snapshots of the fields were taken. From this ensemble, the 
conditional averages have been estimated. Fig. [1] (left) shows the conditional 



averages appearing in (4. 13a I. As can be seen from the figure, the conditional 



average related to the forcing term is approximately linear and thus similar 



to the naive closure assumption (4.19b I. The conditional averages related to 



the diffusive term and the pressure gradient, however, differ strongly from the 
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Figure 5: Comparison of the directly estimated PDF and the one reconstructed using relation 
( |4,15| . The directly estimated PDF is marked by black dots to indicate its presence. 



naive closure. The conditional diffusive term is negatively correlated with the 
velocity, indicating that this term tends to decelerate a fluid particle on average. 
The pressure term shows an interesting zero-crossing, which means that for low 
velocities a fluid particle is on average decelerated whereas for high velocities it 
is accelerated due to this term. For a non-vanishing pressure contribution, the 
occurrence of a zero-crossing is necessary to fulfill the integral constraint (4.161. 
The terms of the conditional dissipation tensor are depicted in fig. [4] (right). It 
becomes clear that the eigenvalues of the conditional dissipation tensor approxi- 
mately coincide, i. e. the conditional dissipation tensor is approximately isotropic. 
For low values of velocity the naive closure assumption (4.201 is acceptable; 
strong deviations, however, occur for high values of velocity indicating that high 
velocity regions in the flow are on average related to regions of increased kinetic 
energy dissipation. As a result, a simple decoupling argument of large-scale and 
small-scale fields yields only a rough approximation to the observed DNS data. 

Fig. [5] shows the evaluation of the stationary solution (4.151 with the con- 
ditional averages estimated from the DNS data. As expected, the agreement 
between the directly estimated PDF and the result of (4.151 is perfect, but 
both PDFs deviate notably from the Maxwellian shape expected for a Gaussian 
distributed velocity vector. Thus the conclusion can be drawn that the observed 
deviations from Gaussianity of the single-point velocity PDF are the result of a 
rather intricate interplay of statistical correlations. 



5. Vorticity Statistics in Three-Dimensional Turbulence 

5.1. LMN Hierarchy for the Vorticity 

Analogous to the LMN hierarchy for the velocity one can derive a hierarchy 
of evolution equations for the vorticity PDFs, which first has been formulated by 
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Novikov |12| . The hierarchy is based on the three-dimensional vorticity equation 



d 

+ u(x,t) ■ V )u)(x,t) = S(x,t)u(x,t) + uAu(x,t) + V x F(x,t) (5.1) 



dt 



where S = ^[Vu + (Vm) ] is the rate-of-strain tensor. The velocity field can be 
expressed as 

u(x,t) = fdx'u(x',t) x K(x - x') (5.2) 

using the kernel (4.3 1. In the vorticity equation the velocity and the rate-of-strain 
tensor represent the nonlocal character of the underlying dynamics whereas in 
the Navier-Stokes equation it is represented by the pressure gradient. 

In the same manner as before the evolution equation for the iV-point PDF 
can be derived from the dynamic equation for the corresponding quantity. In 
the case of the vorticity we obtain 

8 N f 

— f{{u> h xi},t) + Y/J X] ■ jdx'du'u' x K{xj -x')f(u',x',{ui,xi},t) = 

N 

-^V Wj . • / dx'du>' [J x K(xj -x')]uj ■ V x >f(u)',x',{u>t,xt},t) 

3=1 
N 

- ' / dx ' duJ ' LJ ' S ( x J - x')A x ,f{J, x', {wj, xi}, t) . (5.3) 

3=1 

The first term on the right-hand side stems from the vortex stretching term. 
The forcing term was omitted for the sake of simplicity. Again we see a coupling 
of the iV-point PDF equation to the (N + l)-point equation and obtain therefore 
an infinite hierarchy of coupled evolution equations. 

It is interesting to note here that the term resulting from the rate-of-strain 
tensor in eq. (5.3 1 is related to the conditional first-order moment of vorticity. 
This has to be contrasted with the fact that the term due to the pressure gradient 
in eq. (4.4 1 is related to a conditional second-order moment. Therefore, loosely 
speaking, less information is needed to close the PDF equation for the vorticity 
than the one for the velocity. The conditional vorticity field has been investigated 
in [721173]. 

An alternative representation of the LMN hierarchy for the vorticity is 
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obtained when we express the terms in (5.3 1 as conditional averages: 



a N 

~f({ut,xi},t) + • (u{xj,t)\{u} h xi})f({uji,xi},t) 

3 = 1 
N 



3 = 1 
N 

3 = 1 
N 

-X] V «i-( V ^ xFfo.i)^,,^})/^,,^},*) (5.4) 

3 = 1 

In this way we again obtain a formally closed equation, which however contains 
the unknown conditional averages. These averages can be either modeled or 
estimated from experiments or simulations. 



Eq. (5.4 1 shows that the statistics of the rate-of-strain tensor, the symmetric 
part of the velocity gradient tensor, plays an important role for the statistics 
of the vorticity, the anti-symmetric part. Thus it is not surprising that the 
statistics of the full velocity gradient tensor has received a lot of interest lately. 
We refrain here from going into the details and refer the interested reader to |74"] 
for an overview over a number of recent modeling approaches and to [T] for a 
discussion of various statistical properties of the velocity gradient tensor. 

5.2. Single-Point Vorticity Statistics 

In section 14.21 we have demonstrated how a combination of the kinetic 
equation with DNS results has led to insights regarding the slight deviations 
from Gaussianity of the single-point velocity PDF. Obviously this analysis is 
also feasible for the vorticity PDF. However, the results are expected to differ 
strongly from the velocity case as the vorticity PDF is known to be far from 
Gaussian with an enhanced probability to encounter strong vorticities in the 
flow. This observation can be seen as the statistical reflection of the occurrence 
of coherent vortex structures in the flow (see fig. [3| that leads to pronounced 
spatio-temporal correlations. 

The kinetic equations for the single-point velocity PDF have been discussed 
comparably detailed above, and the derivation is completely analogous for the 
vorticity. Hence we focus stronger on a discussion of the DNS results and a 
comparison of the velocity and vorticity statistics in this section. We kindly refer 
the reader to [7SJ H3] for a detailed presentation. 

Introducing the conditional averages on the single-point level leads to a 



kinetic equation completely analogous to (4.6 1: 
d r i r 

—f(u,x,t)+V- |(u|u;)/(u;,a:,i)J = -V w • |^(Su; + ^Au; + V x F\u)f{u),x,t) 

(5.5) 
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Here appear the advective derivative of the PDF on the left-hand side and 
the conditionally averaged terms of the vorticity equation on the right-hand 
side, implying that the evolution of the PDF is governed by the statistics of 
vortex stretching, vorticity diffusion and the external forcing. If we assume a 
homogeneous flow, the advective term on the left-hand side vanishes as neither 
the conditional velocity nor the PDF itself depend on the spatial coordinate. 

As in the velocity case, we also exploit statistical isotropy. While the diffusive 
term and the forcing term are treated like the corresponding terms from the 
Navier-Stokes equation, the conditional vortex stretching term requires special 
scrutiny. As we condition on the single-point vorticity, the vorticity realization 
can be replaced by the sample-space vorticity and pulled out of the average, 

(Sw|w> = (S|w)w, (5.6) 

such that the conditional rate-of-strain tensor has to be specified. Here we have 
to take into account that this tensor is traceless. This leads to the following 
general functional form of the unclosed terms: 

(S ij \u) = ^E(u)(3'^--S ij y j = (<2i ■ Suj\u>) (5.7a) 

(z/Au>|u>) = A(cj) Q A(w) = (yd ■ Au>|u>) (5.7b) 

(V x F\u>) = $(w) <2 $H = (u ■ (V x F)\u) (5.7c) 

By this we have again reduced the unclosed terms to scalar quantities depending 
on the scalar argument uj. The conditional rate-of-strain tensor is fully specified 
by one of its eigenvalues E(w). 

Due to homogeneity, we may replace the conditional diffusive term with the 
so-called conditional enstrophy dissipation tensor, which is again characterized 
by its eigenvalues, 

A» = (-f^f^) = M(«) % + [AM - ^ ■ (5-8) 
This tensor contains the statistical correlations of the vorticity with the vorticity 



gradients. Inserting these results into the original kinetic equation (5.5 1 lets 
us write down two versions of kinetic equations for the PDF of the vorticity 
magnitude, 

f(u,t)=-^(^u,t)u + A(u,t) + ^u,t))f(u,t) (5.9a) 



dt ' du 

f) = -~ (e(«, t)u + t) - 2 J*^A^ t) - *)/(«, t) . 

(5.9b) 



For a further discussion of the stationary vorticity statistics, we turn to 
the DNS results, which have been obtained from the DNS data set described 
in section |4.2| The conditional averages arising in the kinetic equations are 
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Figure 6: Conditional averages and PDF for the single-point vorticity statistics. The PDF 
directly estimated from data is marked by black dots to indicate its presence as it is exactly 
overlapped by the reconstructed PDF. 



shown in fig. [6] In the left graph, the conditional vortex stretching term, the 
diffusive term and the term related to the external forcing are shown. The 
vortex stretching term exhibits a nonlinear functional form indicating strong 
statistical correlations of the rate-of-strain-tensor and the vorticity. As the rate- 
of-strain field and the vorticity field are expected to vary on comparable scales, 
a simple scale separation argument as discussed for the velocity does not hold. 
As becomes also clear from fig. [6] (left), the conditional diffusive term almost 
perfectly balances the vortex stretching term. Considering the fact derived from 



the kinetic equation (5.9a) that the conditional balance 

= S(w)w + A(lj) + $(w) (5.10) 

has to hold then leads to the conclusion that the conditional forcing term has 
to vanish approximately, which is numerically confirmed in fig. [6] (left). We 
would like to mention that this result has already been theoretically predicted 
by Novikov |73| on the basis of a high Reynolds number argument. As a result, 
the single-point statistics of the vorticity is rather independent of the external 
forcing and the approximate relation 

« S(w) uj + A(uj) (5.11) 

holds. In this sense, this statistical quantity purely reflects the nonlinear vor- 
ticity dynamics, and it is tempting to study the impact of the coherent vortex 
structures on this statistical balance. This has been done with the help of wavelet 
decomposition techniques in ref. [75] to which we kindly refer the interested 
reader. 

Also shown in fig. [6] (right) are the eigenvalues of the conditional enstrophy 
dissipation tensor, which display a strong w-dependence. This is in contrast to 
the conditional energy dissipation tensor discussed in section [472] There, the 
simple scale-separation argument between the velocity and (pseudo-)dissipation 
fields has turned out to be approximately valid for low velocities. For the 
vorticity, much stronger correlations are observed because the vorticity and the 
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gradients of the vorticity vary on comparable scales. Also the two eigenvalues 
differ strongly in amplitude (note that fj, and A/4 are shown in fig. [6] (right)). 
This is an indication of strong directional correlations of the vorticity and the 
vorticity gradients related to the filamentary coherent vortex structures. It is 



again possible to calculate the stationary solution of (5.9b) which takes the form 



>< w > = w) exp L duJ m ■ (5 - 12) 

As we have learned from the conditional balance of vortex stretching and vorticity 
diffusion, the forcing term $ can to good approximation be neglected in this 
relation. This is a quantitative proof of the above statement that the single-point 
vorticity statistics is rather independent of the external forcing mechanism. As 
pointed out by Novikov [73 , the influence of the forcing is rather implicit as 
it maintains the stationarity of the flow. Novikov, however, also showed that 
the external forcing is expected to have influence on the two-point vorticity 



correlations [77] . The validity of relation (5.121 is again tested with our numerical 



results and, as expected, excellent agreement is demonstrated in fig. [6j 

To summarize the results of the investigation of the three-dimensional single- 
point velocity and vorticity statistics, by the joint use of kinetic equations 
and DNS data it has become possible to account the slight deviations from 
Gaussianity of the single-point velocity to moderate statistical correlations and 
the strong non-Gaussianity of the vorticity PDF to more pronounced statistical 
correlations. The fact that the vorticity statistics is more independent of the 
large-scale forcing together with the occurrence of coherent vortex structures 
indicates that a deeper understanding of the statistics may result from further 
investigation of vortex dynamics. 



6. Application to Turbulent Rayleigh-Benard Convection 

Kinetic equations have been investigated for Rayleigh-Benard convection by 
Lulff et al. [78 . Already the single-point kinetic equation yields considerable 
insight into the connection between statistics and dynamics of turbulent con- 
vection, and also allows one to link to phenomenological theories formulated 
on the basis of dimensional arguments. Especially, basic assumptions of the 
Grossmann-Lohse theory [75] can be investigated by the combination of kinetic 
equations closed by taking conditional averages from DNS. 

Rayleigh-Benard convection idealizes the heat transport of a convecting fluid 
enclosed between a hot bottom and a cold top plate. It is described by the 
Oberbeck-Boussinesq equations for the velocity field u(x, t) and the temperature 
field T(x,t): 

d 

—T + u ■ VT = AT 
at 

d 

— u + u- Vit = -Vp + Pr Aw + PrRaTe^ (6.1) 
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Figure 7: Visualization of the temperature field in turbulent Rayleigh-Benard convection. Red 
represents hot, blue cold parts of the fluid. Temperatures around the mean are translucent. 
Parameters are Ra = 2.4 X 10 7 , Pr = 1, aspect ratio T = 4. Videos of this visualization to be 
found at 1621 . 



The velocity field is incompressible, V • u — 0, and vanishes at solid boundaries. 
Additionally, boundary conditions for the temperature field have to be defined, 
i. e. a hot bottom and a cold top plate. Gravity acts in vertical e 2 -direction. The 
system is characterized by two dimensionless parameters, the Rayleigh number 
Ra, which describes the temperature gradient, and the Prandtl number Pr, which 
describes fluid properties. 

To estimate conditional averages from DNS, the full Oberbeck-Boussinesq 
equations are solved via a pseudospectral method acting on an equidistant, 
tri-periodic grid. The boundary conditions in vertical direction, i. e. impermeable 
plates of constant temperature, are enforced via a suitably designed volume 
penalization approach (for more information about the method of volume penal- 
ization, we refer the reader to e. g. jSHl EU E2 ) , while in the horizontal directions, 
periodic boundary conditions are assumed. A volume rendering of the turbulent 
temperature field is exhibited in fig. [7J displaying the overall tangled structure 
of the temperature field, as well as thin filaments of hot and cold fluid emerging 
from the horizontal plates which are referred to as plumes. 



6.1. LMN Hierarchy for Temperature and Velocity 

The full single-point kinetic equation for the joint PDF of temperature and 
velocity, f(T,u,x,t), has been derived in |78| : 



d d 



AT\T,u,x,t)f 
((-Wp + PrAu\T,u,x,t) + PrRaTe z )/ 



■2) 
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For Pr = 1, a homogeneity relation similar to (4.7 1 can be applied, and the 
resulting kinetic equation reads 



d 
Of 



f + u-Vf 



A/-V, 



dTdu 



«-Vp|*)+RaTe,)/ 
(VT-Vuj-I*)/ 



2 r 



d 2 

duidu 



((VT) 2 K)/ 



(6.3) 



with * = T,u,x,t abbreviating the dependencies of the conditional averages. The 
joint temperature-velocity PDF is effectively determined by the conditionally 
averaged dissipation-like terms ((VT) 2 |*), (VT • Vuj\*) and (Vitj • Vitj|*), and 
the conditional pressure gradient (— Vp|*). 

From the full kinetic equation (6.2 1 one can proceed to the reduced PDF 



h(T, z) of temperature alone, which after applying the statistical symmetries of 
homogeneity and stationarity depends on the vertical coordinate z only and is 
governed by the kinetic equation 



d_ 

dz 



(u z \T, z)h{T,z) 



d_ 
df 



(AT\T, z)h(T,z) 



= 0. 



(6.4) 



Thereby, the average dynamics in T, z-phase space and hence the temperature 
PDF is determined by the conditional averages of vertical velocity (u z \T, z) and 
heat diffusion (AT\T, z). The associated characteristic equations that depict the 
aforementioned dynamics read 



d 
d^ 

_d_ 

d^ 

d 



T(s) = (AT\T,z) 
z(s) = (u z \T,z) 



T=T(s) 
z—z{s) 

T=T(s) 
z—z(s) 



h(s) 



d 



'u~\T, z) 



d 



(AT\T,z) 



t=t(s) M s ) ■ 

z=z(s) 



(6.5a) 
(6.5b) 
(6.5c) 



ds w \dz" 21 ' ' dT 
with s being the parametrization of the characteristic curves; s can be identified 



with the time t of the system (cf. section 2.4 and eq. (2.42a I therein) 



6.2. Numerical Results 

In fig. [8] the z-resolved temperature PDF h(T,z) is shown. In the left 
panel, the full two-dimensional PDF is color-coded, and one can observe that 
the PDF shows an almost constant mean value in the bulk regions (i. e. for 
4<5t ;$ z < 1 — A8t) while the skewness varies with z (indeed, it is observed 
that the skewness varies almost linear in the bulk, cf. |83J). In the boundary 
regions the PDF contracts to a (5-function due to the boundary conditions of 
fixed temperature. The right panel of fig. [8] shows one-dimensional cuts at fixed 
values of z as indicated in the two-dimensional PDF on the left; one can see 
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Figure 8: Height- resolved temperature PDF h{T, z) for turbulent Rayleigh-Benard convection 
in logarithmic scaling. Dimensionless parameters as in fig. Left: Two-dimensional color 
plot. The solid black and gray lines indicate mean and standard deviation of the temperature 
PDF, respectively. Right: Slices in T-direction, as indicated by the horizontal dashed lines 
in the color plot on the left, with <5y = being the thermal boundary layer thickness and 

T rms = v (T 2 )v the quadratic mean of temperature with respect to the whole volume. 



that the temperature PDF broadens for low values of z before contracting to a 
S- function. 

The vector field that governs the characteristic curves is depicted in fig. [9] 
together with a color plot of the PDF h(T, z). The characteristic curves show 
the average behavior of a fluid parcel, i. e. the typical RB cycle of fluid heating 
up at the bottom, rising up to the top plate, cooling down at the top and falling 
down to the bottom again. It is tempting to interpret the characteristics as a 
kind of Lagrangian dynamics of a tracer particle inside the RB cell, in close 
analogy to the smart particles experimentally investigated by the Lyon group 
(see e. g. [HI])- However, in contrast to Lagrangian trajectories the characteristics 
are trajectories in an averaged field and, therefore, only take the chaotic features 
into account in an averaged sense. 

The alternative representation of the heat conductive term similar to the 



homogeneity relation which leads to (6.3 1 gives rise to the kinetic equation 



d_ 

dz~ 



{u z \T, z)h(T,z] 



a 2 
dz 



2 h(T,z) 



d 2 

Q T 2 



((VT) 2 \T,z)h(T,z) 



(6.6) 



Here, the conditional average ((Vrr|T, z) conies up. This term can be related 
to the Nusselt number, which is defined as the spatial average over the whole 
cell with volume V, Nu = ((VT) 2 )^. This quantity can be calculated in terms 



of the conditional average z) as 

yfdx fdT((VT) 2 \T,z) h(T, z). 



Nu 



(6.7) 



Therefore, ((VTy\T,z) can be viewed as a conditional Nusselt number. It is of 
considerable interest for the evaluation of phenomenological theories concerning 
the Rayleigh number dependence of the Nusselt number based on a decomposition 
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Figure 9: Vector field (T(s), z(s)) governing the characteristic equations in turbulent Rayleigh- 
Benard convection, cf. eq. l |6.5[ l. Additionally, the temperature PDF is color-coded in compliance 
with fig. [8] Dimensionless parameters as in fig. [7] One can clearly see the generic RB cycle 
discussed in the text. 



of the heat transport into bulk and boundary contributions, which underlies 
the Grossmann-Lohse theory |79| . Also, this theory relies on the z-resolved 
dissipation rates, which can be linked to the conditional averages of quadratic 



terms in (6.3 I which naturally come up in our derivation. 



An approach similar to the one pursued by us has been performed by Yakhot 
|18) and Ching [25] for the temperature PDF, and also by Sinai et al. |17| 
for a passive scalar (in contrast to the temperature as an active scalar in RB 
convection). Starting from different premises than we do, Yakhot ends up with a 
formula that also describes the temperature PDF in terms of conditional averages. 
After approximating these averages in the high Rayleigh number limit, he arrives 
at a temperature PDF of exponential shape in the whole fluid volume so that 
the PDF does not depend on the vertical coordinate z. In contrast to this result 
and also to summarize our findings, the temperature PDF that we measured 
from the numerics clearly shows a z-dependence, while also our ansatz is able 
to connect the statistics of the system, i. e. the different conditional averages 
that naturally show up in our derivation, to the dynamics of the system, i.e. the 
characteristic curves that reproduce the distinctive RB cycle. 



7. Burgers Turbulence 

Up to now we have reviewed kinetic equations, which were closed using con- 
ditional expectations taken from DNS. In this section, we review an application 
of the LMN hierarchy to the forced Burgers equation, which allows for an exact 
assessment of the multi-point PDFs, and hence yields considerable insights into 
the structure of the LMN hierarchy. Burgers turbulence has always been used 
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as a playground for the evaluation of novel theoretical approaches and concepts 
to the problem of turbulence [29 . 

The Burgers equation can be exactly solved by application of a Hopf-Cole 
transformation. However, the statistical behavior of the forced Burgers equation 
is still under heavy debate. We refer the reader to the review of Bee and Khanin 



Polyakov [25] investigated the increment statistics of the driven Burgers 
equation 

d 2 

u(x,t) = v—ru(x,t) + F(x,t) (7.1) 



d . . d 
— + u{x,t) — 
at ox 



dx 2 



with a long-range correlated white noise source, based on the cumulative 
probability (Q(v — [u(x + r,t) — u(x,i)])) (S(«) denotes the Heaviside function, 
Q'(v) — S(v)). He introduced a closure for the expectation of the conditional 
dissipation term and was able to suggest a PDF for the velocity gradient. This 
work has stimulated subsequent investigations on the treatment of the dissipation 
anomaly. 

In [55] it was possible to find an exact solution of the forced Burgers equation 
with a stochastic forcing F(x, t) = rj(t)x that is linear in x and 77(f) is a Gaussian 
white noise term with (r}(t)rj(t')) — QS(t — t'). The solution has the form 

u(x, t) = a(t)x + A(t)w(A(t)x, r(f)) . (7.2) 

Here, w(£,t) is an arbitrary solution of the homogeneous Burgers equation, and 
the amplitude a(f) obeys the stochastic differential equation a = —a 2 + 77. The 
observables A(f) and r(f) are random variables related to the statistics of the 
amplitude a(t) by A = — aA and r = A 2 . The solution exhibits a finite-time 
singularity, which corresponds to the formation of a new shock. It is interesting 
to notice that the iV-point PDFs can be constructed explicitly in terms of the 
statistics of a(t) and the initial condition for the solution of the homogeneous 
Burgers equation w(£,r). For the sake of convenience, we formulate the TV-point 
PDF as 

f(Ul,X!,. . . , u N , x N ,t) 

r N 

I dadAdr h(a, A,T,t)Yl(5(ui - ax t - Aw(Ax l , t))) w , (7.3) 

^ 4 = 1 

where h(a, A, r, f) is the solution of the Fokker-Planck equation for the stochastic 
process 

a(t) = -a 2 {t)+r](t) 
A(f) = e -Jo td *XO 

r(f) - f dt" e- 2 ti" dt ' Q ( t ' ) . (7.4) 



Here, an additional average (•)„, with respect to the initial conditions of the 
solution w(£,t) has been included. 
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It is interesting to note that, if we take the initial condition w = 0, the 
TV-point PDF has the form 



f(ui,Xx, . . . ,u N ,x N ,t) = 

Slut - — U2)6\U2 - — u 3 ) . . .8\u N -i - XN ~ 1 u N ) — h (—,t] , 
V x N J \ x N J \ x N J x N \x N J 

(7.5) 

where ho is the solution of the Fokker-Planck equation 
d d Q d 2 

j-h (a,t) = —a 2 h (a,t) + ——-h (a,t) . (7.6) 

Although the result is rather trivial, due to the trivial behavior of u(x,t), we 
can conclude that /(ui, X\, . . . , upf, xjq, t) actually describes a Markov chain in 
scale with a sharp transition probability 

p(ui,Xi\u i+ i,Xi + i) =S[Ui - X ^u i+1 ) . (7.7) 
V x N J 

We remind the reader that such a structure underlies the phenomenological 
theory of intermittency formulated by Friedrich and Peinke [30] for the direct 
cascade of turbulence. It would be interesting to see whether the Markovian 
property, which can be proven here from first principles, survives for the case, 
where the forcing F(x, t) has a finite correlation length leading to the emergence 
of shocks. 



8. Conclusions 

We have presented an overview of recent works on the framework of the 
Lundgren-Monin-Novikov hierarchy. After giving a short introduction into the 
basic methods and concepts, we have reviewed applications to various turbulent 
systems ranging from the vorticity statistics in two-dimensional turbulence over 
the single-point velocity and vorticity statistics in three-dimensional turbulence 
to the temperature statistics in turbulent Rayleigh-Benard convection. In 
these works, the closure problem has been treated by estimating the unclosed 
terms from direct numerical simulations. While this approach does not solve 
the closure problem, a direct link between basic dynamical features of the 
system under consideration and the observed statistics is established without 
further phenomenological assumptions. By this, insights into various hallmarks 
of turbulence statistics like the origin of the energy and enstrophy transfer 
across scales in two-dimensional turbulence, deviations from Gaussianity in 
three-dimensional turbulence and the height-resolved temperature statistics in 
turbulent Rayleigh-Benard convection have been obtained. Furthermore, for the 
linearly random driven Burgers equation we have reviewed an exact solution for 
the LMN hierarchy. 
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An interesting point is the fact that the characteristic equations associated 
to the kinetic equation for the A-point Eulerian PDF can be interpreted in a 
Lagrangian sense as the description of particles moving in an averaged field. 
Especially in two-dimensional turbulence, this brings up the notion of screened or 
quasi- vortices, which contain the statistically averaged structure of the vorticity 
field. In the notion of fluid dynamics this directly links the LMN hierarchy 
to the field of subgrid modeling. Furthermore, it allows to make contact with 
Lagrangian models for turbulence like, for example, the tetrad models, reviewed 
by Pumir and Naso in this volume. 

In the future, the investigation of the PDF of velocity increments started by 
Yakhot |19| should be complemented by an estimation of the arising conditional 
expectations of pressure gradient and dissipation from DNS. An assessment 
of the corresponding kinetic equation would allow to make contact with the 
phenomenological theories of the turbulent cascade based on the multifractal 
idea (see for example the review of Chevillard et al. in this volume) or the 
phenomenological theory of Friedrich and Peinke [30 , based on the stochastic 
analysis of the statistics of velocity increments. 

A shortcoming of the LMN approach clearly originates from the fact that an 
estimation of conditional PDFs from DNS has up to now only been feasible for 
one- or two-point PDFs. It would be interesting to investigate whether a suitable 
application of more sophisticated methods of data analysis could help to extend 
the approach to more points leading to a formulation of higher-order kinetic 
equations. Another possibility to overcome this restriction is the development 
of approximate analytical expressions for the conditional expectations. Up to 
now, conditional expectations based on Gaussian statistics have been used, 
which explicitly fail for the two-dimensional inverse cascade as well as the single- 
point PDF of three-dimensional vorticity. However, the detailed investigation 
of the failure sheds considerable light on the necessary extension of the simple 
conditional Gaussian ansatz. Both extensions are expected to help to advance 
significantly our understanding of the statistics of fully developed turbulence. 
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